Microscopic statistical basis of classical thermodynamics of finite systems. 
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Heat can flow from cold to hot at any phase separation. Therefore Lynden-Bell's gravo-thermal 
catastrophe [lj must be reconsidered. The original objects of Thermodynamics, the separation 
of phases at first order phase transitions, like boiling water in steam engines, are not described 
by a single canonical ensemble. Inter-phase fluctuations are not covered. The basic principles of 
statistical mechanics, especially of phase transitions have to be reconsidered without the use of the 
thermodynamic limit. Then thermo-statistics applies also to nuclei and large astronomical systems. 
A lot of similarity exists between the accessible phase space of fragmenting nuclei and inhomogeneous 
multi stellar systems. 
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Since the beginning of Thermodynamics in the first 
half of the 19. century its original motivation was the de- 
scription of steam engines and the liquid to gas transition 
of water. Here water becomes inhomogeneous and devel- 
ops a separation of the gas phase from the liquid, i.e. 
water boils. This will be analyzed from the new perspec- 
tive of microcanonical statistics in section ITTT1 

A little later statistical mechanics was developed by 
BoltzmannQ to explain the microscopic mechanical ba- 
sis of Thermodynamics. Up to now it is generally be- 
lieved that this is given by the Boltzmann-Gibbs canon- 
ical statistics. As traditional canonical statistics works 
only for homogeneous, infinite systems, phase separations 
remain outside of standard Boltzmann-Gibbs thermo- 
statistics, which, consequently, signal phase-transitions 
of first order by Yang-Lee singularities. 

It is amusing that this fact that is essential for the orig- 
inal purpose of Thermodynamics to describe steam en- 
gines was never treated completely in the past 150 years. 
The system must be somewhat artificially split into (still 
macroscopic and homogeneous) pieces of each individ- 
ual phase Q . The most interesting configurations of two 
coexisting phases cannot be described by a single canoni- 
cal ensemble. Important inter-phase fluctuations remain 
outside, etc. This is all hidden due to the restriction to 
homogeneous systems in the thermodynamic limit. 

Also the second law can rigorously be formulated only 
microcanonically: Already Clausius |4| distinguished be- 
tween external and internal entropy generating mecha- 
nisms. The second law is only related to the latter mech- 
anism |{| , the internal entropy generation. Again, canon- 
ical Boltzmann-Gibbs statistics is insensitive to this im- 
portant difference. 

For this purpose, and also to describe small systems 
like fragmenting nuclei or non-extensive ones like macro- 
scopic systems at phase-separation, or even very large, 
self-gravitating, systems, we need a new and deeper def- 
inition of statistical mechanics and as the heart of it: of 
entropy. For this purpose it is crucial to avoid the ther- 



modynamic limit. 



I. WHAT IS ENTROPY? 

Entropy, S, is the characteristic entity of thermody- 
namics. Its use distinguishes thermodynamics from all 
other physics; therefore, its proper understanding is es- 
sential. The understanding of entropy is sometimes ob- 
scured by frequent use of the Boltzmann-Gibbs canon- 
ical ensemble, and the thermodynamic limit. Also its 
relationship to the second law is beset with confusion be- 
tween external transfers of entropy d e S and its internal 
production diS. 

The main source of the confusion is of course the lack 
of a clear microscopic and mechanical understanding of 
the fundamental quantities of thermodynamics like heat, 
external vs. internal work, temperature, and last not 
least entropy, at the times of Clausius and possibly even 
today. 

Clausius defined a quantity which he first called 
the "value of metamorphosis" , in German " Verwand- 
lungswert" in 4f]. Eleven years later he [|| gave it the 
name "entropy" S: 



dE 

-jT, 



(1) 



where T is the absolute temperature of the body when 
the momentary change is done, and dE is the increment 
(positive resp. negative) of all different forms of energy 
(heat and potential) put into resp. taken out of the sys- 
tem. 

From the observation that heat does not flow from cold 
to hot (see section [Hi formula however section UTTIl he 
went on to enunciate the second law as: 



AS" 



dE 



>0, 



(2) 
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which Clausius called the uncompensated metamorphosis. 
As will be worked out in section II I II the second law as 
presented by eq.(|2J) remains valid even in cases where 
heat flows from low to higher temperatures. 
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Prigogine c.f. 0, quite clearly stated that the vari- 
ation of S with time is determined by two, crucially dif- 
ferent, mechanisms of its changes: the flow of entropy 
d e S to or from the system under consideration; and its 
internal production diS. While the first type of entropy 
change d e S (that effected by exchange of heat d e Q with 
its surroundings) can be positive, negative or zero, the 
second type of entropy change diS is fundamentally re- 
lated to its spontaneous internal evolution ("Verwand- 
lungen", "metamorphosis" 0) of the system, and states 
the universal irreversibility of spontaneous transitions. It 
can be only positive in any spontaneous transformation. 

Clausius gives an illuminating example in When 
an ideal gas suddenly streams under isolating conditions 
from a small vessel with volume V\ into a larger one 
(Va > Vx), neither its internal energy U, nor its tem- 
perature changes, nor external work done, but its in- 
ternal (Boltzmann-)entropy Si eq.© rises, by AS = 
iVln(l^/Vi) . Only by compressing the gas (e.g. isen- 
tropically) and creating heat AE = £a[(F 2 /Vi) 2/3 - 1] 
(which must be finally drained) it can be brought back 
into its initial state. Then, however, the entropy change 
in the cycle, as expressed by integral @, is positive 
(= TV In (V2/V1)). This is also a clear example for a mi- 
crocanonical situation where the entropy change by an 
irreversible metamorphosis of the system is absolutely 
internal. It occurs during the first part of the cycle, the 
expansion, where there is no heat exchange with the en- 
vironment, and consequently no contribution to the inte- 
gral©- The construction by eq.|0) is correct though arti- 
ficial. After completing the cycle the Boltzmann-entropy 
of the gas is of course the same as initially. All this 
will become much more clear by Boltzmann's microscopic 
definition of entropy, which will moreover clarify its real 
statistical nature: 

BoltzmannQ later defined the entropy of an isolated 
system (for which the energy exchange with the environ- 
ment d e Q = 0) in terms of the sum of possible configura- 
tions, W, which the system can assume consistent with 
its constraints of given energy and volume: 



S=k*lnW 



(3) 



as written on Boltzmann's tomb-stone, with 



W(E,N, V) 



d 3N p d 3N q 

N\(2irh) 3N 



e 6{E-H{q,P}) 



in semi-classical approximation. E is the total energy, N 
is the number of particles and V the volume. Or, more 
appropriate for a finite quantum-mechanical system: 



W(E, N,V) 



Eall cigenstates n of H with given N,V, /r\ 
and E < E„ < E + e l°J 



and eo ~ the macroscopic energy resolution. This is still 
up to day the deepest, most fundamental, and most sim- 
ple definition of entropy. There is no need of the ther- 
modynamic limit, no need of concavity, extensivity and 



homogeneity. In its semi-classical approximation, eq.Q, 
W(E, N, V, ■ ■ ■ ) simply measures the area of the sub- 
manifold of points in the 6 TV-dimensional phase-space (r- 
space) with prescribed energy E, particle number TV, vol- 
ume V, and some other time invariant constraints which 
are here suppressed for simplicity. Because it was Planck 
who coined it in this mathematical form, I will call it the 
Boltzmann-Planck principle. It is further important to 
notice that S(E, TV, V) is everywhere analytical in E 0. 
In the microcanonical ensemble are no "jumps" or multi- 
valuedness in S(E), independently of whether there are 
phase transitions or not, in clear contrast to the canoni- 
cal S(T, N,V). A fact which underlines the fundamental 
role of microcanonical statistics. 



The Boltzmann-Planck formula has a simple but deep 
physical interpretation: W or S measure our ignorance 
about the complete set of initial values for all 6TV mi- 
croscopic degrees of freedom which are needed to spec- 
ify the TV-body system unambiguously 8] . To have com- 
plete knowledge of the system we would need to know 
(within its semiclassical approximation Q) the initial 
positions and velocities of all TV particles in the system, 
which means we would need to know a total of 6TV values. 
Then W would be equal to one and the entropy, S, would 
be zero. However, we usually only know the value of a 
few parameters that change slowly with time, such as the 
energy, number of particles, volume and so on. We gener- 
ally know very little about the positions and velocities of 
the particles. The manifold of all these points in the 6TV- 
dim. phase space, consistent with the given macroscopic 
constraints of E, TV, V, ■ ■ ■ , is the microcanonical ensem- 
ble, which has a well-defined geometrical size W and, by 
equation ©, a non- vanishing entropy, S(E, TV, V, ■ ■ ■ ). 
The dependence of S(E, TV, V, • • • ) on its arguments de- 
termines completely thermostatics and equilibrium ther- 
modynamics. 



Clearly, Hamiltonian (Liouvillean) dynamics of the 
system cannot create the missing information about the 
initial values - i.e. the entropy S(E, TV, V, ■ ■ ■ ) cannot de- 
crease. As has been further worked out in 0, the 
inherent finite resolution of the macroscopic description 
implies an increase of W or S with time when an external 
constraint is relaxed. Such is a statement of the second 
law of thermodynamics, which requires that the internal 
production of entropy be positive for every spontaneous 
process. Analysis of the consequences of the second law 
by the microcanonical ensemble is appropriate because, 
in an isolated system (which is the one relevant for the 
microcanonical ensemble), the changes in total entropy 
must represent the internal production of entropy, see 
above, and there are no additional uncontrolled fluctuat- 
ing energy exchanges with the environment. 
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II. THE ZERO'TH LAW IN CONVENTIONAL 
EXTENSIVE THERMODYNAMICS 

In conventional (extensive) thermodynamics thermal 
equilibrium of two systems (1 & 2) is established by 
bringing them into thermal contact which allows free en- 
ergy exchange. Equilibrium is established when the total 
entropy 

S 1+2 (E,E 1 ) = S 1 {E 1 )+S 2 (E-E 1 ) (6) 

is maximal. Under an energy flux AE 2 ^\ from 2 — > 1 
the total entropy changes to lowest order in AE by 

A5i +2 |b = (T 2 -Ti)A^i. (7) 

Consequently, a maximum of Stotal(E, Ei)\e > £1+2 will 
be approached when 



sign(AS', 



total ) 



sign(T 2 - Ti)sign(A£ 2 ^i) > 0. (8) 



From here Clausius' first formulation of the Second Law 
follows: "Heat always flows from hot to cold". Essential 
for this conclusion is the additivity of S under the split 
(eq|SJ). There are no correlations, which are destroyed 
when an extensive system is split. Temperature is an 
appropriate control parameter for extensive systems. 

It is further easy to see that the heat capacity of an 
extensive system with S(E,N) — 2S(E/2, N/2) is neces- 
sarily positive 



Cv(E) = dE/dT 



{dS/OEf 
d 2 S/dE 2 



> 



(9) 



The combination two pieces of N/2 particles each, but 
with different energy per particle, one at e a — e 2 — Ae/2 
and a second at e b — e 2 + Ae/2, must lead to S(E 2 , N) > 
S{E a /2,N/2) + S(E b /2,N/2), the simple algebraic sum 
of the individual entropies because by combining the two 
pieces one normally looses information. This, however, 
is equal to [S(E a ,N) + S(E b ,N)}/2, thus S(E 2 ,N) > 
[S(E a ,N) + S{E b ,N)}/2. I.e. the entropy S(E,N) of 
an extensive system is necessarily concave, d 2 S/dE 2 < 
and eq. ED follows. In the next section we will see that 
therefore extensive systems cannot have phase separation, 
the characteristic signal of transition of first order. 



III. NO PHASE SEPARATION WITHOUT A 
CONVEX, NON-EXTENSIVE S(E) 

At phase separation the weight e s ( E )— E / T of the con- 
figurations with energy E in the definition of the canon- 
ical partition sum 



Z(T) 



e S(E)-E/T dE 



(10) 



becomes bimodal, at the transition temperature it has 
two peaks, the liquid and the gas configurations which 



are separated in energy by the latent heat. Consequently 
S{E) must be convex (d 2 S/dE 2 > 0, like y — x 2 ) and 
the weight in 1)10(1 has a minimum at E m i n between the 
two pure phases. Of course, the minimum can only be 
seen in the microcanonical ensemble where the energy 
is controlled and its fluctuations forbidden. Otherwise, 
the system would fluctuate between the two pure phases 
by an, for macroscopic systems even macroscopic, energy 
AE ~ Ei a t of the order of the latent heat. The heat 
capacity is 



C v {Emin) = dE/dT : 



{dS/dEf 
d 2 S/dE 2 



< 0. 



(11) 



I.e. the convexity of S(E) and the negative heat ca- 
pacity are the generic and necessary signals of phase- 
sevarationflftf. It is amusing that this fact, which is es- 
sential for the original purpose of Thermodynamics to de- 
scribe steam engines, seems never been really recognized 
in the past 150 years. However, such macroscopic energy 
fluctuations and the resulting negative specific heat are 
already early discussed in high-energy physics by Carlitz 

The existence of the negative heat capacity at phase 
separation has a surprising but fundamental conse- 
quence: Combining two equal systems with negative heat 
capacity and different energy per particle, they will re- 
lax with a flow of energy from the lower to the higher 
temperature! This is consistent with the naive picture 
of an energy equilibration. Thus Clausius' "energy flows 
always from hot to cold", i.e. the dominant control-role 
of the temperature in thermo- statistics as emphasized by 
Hertz ]1S] is violated. Of course this shows quite clearly 
that unlike to extensive thermodynamics the temperature 
is not the appropriate control parameter in non-extensive 
situations like e.g. at phase separations, nuclear frag- 
mentation, or stellar systems. \lAl 

By the same reason the well known paradox of Antonov 
in astro-physics due to the occurrence of negative heat 
capacities must be reconsidered: By using standard ar- 
guments from extensive thermodynamics Lynden-Bcll 
[l| claims that a system a with negative heat capacity 
C a < in gravitational contact with another b with pos- 
itive heat capacity C b > will be unstable: If initially 
T a > T b the hotter system a transfers energy to the colder 
b and by this both become even hotter! If C b > —C a , 
T a rises faster than T b and if the heat capacities don't 
change, this will go for ever. This is Lynden-Bells gravo- 
thermal catastrophe. This is wrong because just the op- 
posite happens, the hotter a may even absorb energy from 
the colder b and both systems come to equ ilibrium at the 
same intermediate temperature c.f. [iH Il5j. Negative 
heat can only occur in the microcanonical ensemble. 

As phase separation exists also in the thermodynamic 
limit, by the same arguments as above the curvature of 
S(E) remains convex, d 2 S/(dE) 2 > 0. Consequently, 
the negative heat capacity should also be seen in ordinary 
macroscopic systems studied in chemistry! 

Searching for example in Guggenheims book 3] one 
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finds some cryptic notes in §3 that the heat capacity of 
steam at saturation is negative. No notice that this is 
the generic effect at any phase separation! Therefore let 
me recapitulate in the next section how chemists treat 
phase separation of macroscopic systems and then point 
out why this does not work in non-extensive systems like 
fragmenting nuclei, at phase separation in normal macro- 
scopic systems, or large astronomical systems. 



The positiveness of any canonical Cy (T) or Cp(T) l|15(l 
is of course the reason why the inhomogeneous system of 
several coexisting phases (ai&a^) with an overall nega- 
tive heat capacity cannot be described by a single com- 
mon canonical distribution |llL . 

This new but fundamental interpretation of thermo- 
statistics was introduced to the chemistry community in 

0G3 



IV. MACROSCOPIC SYSTEMS IN CHEMISTRY 



V. NEW KIND OF PHASES WELL SEEN IN 
HOT NUCLEI OR MULTI-STAR SYSTEMS. 



Systems studied in chemical thermodynamics consist 
of several homogeneous macroscopic phases ai,ot2,--- 
cf.0j. Their mutual equilibrium must be explicitly con- 
structed from outside. 

Each of these phases are assumed to be homogeneous 
and macroscopic (in the "thermodynamic limit" (N a — > 
oo\ Pa=CO nst))- There is no common canonical ensemble 
for the entire system of the coexisting phases. Only 
the canonical ensemble of each phase separately becomes 
equivalent in the limit to its microcanonical counterpart. 

The canonical partition sum of each phase a is defined 
as the Laplace transform of the underl ying microcanoni- 
cal sum of states W(E) a = e s "^ Ella 



e S«(E)-E/T adE _ 



The mean canonical energy is 

< E a {T a ) >=T*d\nZ a (T a )/dT a . 



(12) 



(13) 



In chemical situations proper the assumption of homoge- 
neous macroscopic individual phases is of course accept- 
able. In the thermodynamic limit (N a — ► oo| Pa=const ) of 
a homogeneous phase a, the canonical energy 
<E a (T a )> becomes identical to the microcanonical en- 
ergy E a when the temperature is determined by T~ = 
dS a (E,V a )/dE a . The relative width of the canonical 
energy is 



AE(T) C 



V< El > T - < E a > 



1 



< E a > T 



oc 



The heat capacity at constant volume is 



v a — 



< K >T a ~<E a >\ 
J>2 



> 0. 



(14) 



(15) 



Only in the thermodynamic limit (N a — > oo\ Pci=cons t) 
does the relative energy uncertainty AE a — > 0, and the 
canonical and the microcanonical ensembles for each ho- 
mogeneous phase (a) become equivalent. This equiva- 
lence is the only justification of the canonical ensemble 
controlled by intensive temperature T, or chemical poten- 
tial /j,, or pressure P. I do not know of any microscopic 
foundation of the canonical ensemble and intensive con- 
trol parameters apart from the limit. 



The new lesson to be learned is that if one defines the 
phases by individual peaks j2(| in e s ^ E ^ E ^ T in IjlUI) . then 
there exist also inhomogeneous phases like in fragmented 
nuclei or stellar systems. The general concept of thermo- 
statistics becomes enormously widened. 

Now, certainly neither the phase of the whole multi- 
fragmented nucleus nor the individual fragments them- 
selves can be considered as macroscopic homogeneous 
phases in the sense of chemical thermodynamics (ChTh). 
Consequently, (ChTh) cannot and should not be applied 
to fragmenting nuclei and the microcanonical description 
is ultimately demanded. This becomes explicitly clear by 
the fact that the configurations of a multi-fragmented nu- 
cleus have a negative heat capacity at constant volume 
Cy and also at constant pressure Cp (if at all a pres- 
sure can be associated to nuclear fragmentation 11]). 
Meanwhile a huge amount of experimental evidences of 
negative heat capacities has accumulated: Nuclear frag- 
mentation e.g. [20|. atomic clusters e.g. [2l|. astrophysics 
e -g- |22j |. conventional macroscopic systems at phase sep- 
aration e.g.Q. 

The existence of well defined peaks (i.e. phases as de- 
fined above) in the event distribution of nuclear frag- 
mentation data is demonstrated very nicely in |2^| from 
various points of view. A lot more physics about the 
mechanism of phase transitions can be learned from such 
studies. 



VI. OUTLOOK 

It is a deep and fascinating aspect of nuclear fragmen- 
tation: First, in nuclear fragmentation we can measure 
the whole statistical distribution of the ensemble event 
by event including eventual inter-phase fluctuations. Not 
only their mean values are of physical interest. Statisti- 
cal mechanics can be explored from its first microscopic 
principles in any detail well away from the thermody- 
namic limit. Initiated by our theoretical studies of nu- 
clear fragmentation we found the very general appear- 
ance of a backbending caloric curve T(E) corresponding 
to a negative heat capacity, e.g.:[24| [an early review 
in |23] similar effects were proposed in the melting of 
atomic clusters |26|. Years later its existence in nuclear 
fragmentation was verified experimentally |20j. However, 
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the necessary convexity of the entropy S(E) at any phase 
separation seems to be little known in thermodynamics. 
Clausius' version of the second law "heat always flows 
from hot to cold" is in general violated at any phase sepa- 
ration even in macroscopic systems. Nowadays, the non- 
equivalence of the micro- and the canonical ensemble at 
phase-separations is discussed by many authors, see e.g. 
several relevant papers in . 

In nuclear fragmentation there may be other conserved 
control parameters besides the energy: E.g. in the recent 
paper by Lopez et al. the importance of the total 
angular momentum of the excited nucleus was empha- 
sized. In this paper a bimodality, i.e. phase separation, 
in the mass-asymmetry of the fragments is demonstrated 
controlled by the transferred spin and not by excitation 
energy. This is an interesting, though still theoretical, ex- 
ample of the rich facets of the fragmentation phase tran- 
sition in finite systems which goes beyond the liquid-gas 
transition and does not exist in macroscopic chemistry. 



Angular momentum is also a very crucial control param- 
eter in stellar systems, see below. 

Second, and this may be more important: For the first 
time phase transitions to non- homogeneous phases can be 
studied where these phases are within themselves com- 
posed of several nuclei. This situation is very much analo- 
gous to multi star systems like rotating double stars dur- 
ing intermediate times, when nuclear burning prevents 
their final implosion. The occurrence of negative heat 
capacities is an old well known peculiarity of the thermo- 
statistics of self-gravitating systems pi Also these 
cannot be described by a canonical ensemble. It was 
shown in |lft ITfij how the microcanonical phase space 
of these self-gravitating systems has the realistic config- 
urations which are observed. Of course, the question 
whether these systems really fill uniformly this phase 
space, i.e. whether they are interim equilibrized or not is 
not proven by this observation though it is rather likely. 
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